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Abstract 

We report on the measurement of freeze-out radii for pairs of identical-charge pions measured in 
Pb-Pb collisions at y/sm=2.76 TeV as a function of collision centrality and the average transverse 
momentum of the pair kj. Three-dimensional sizes of the system (femtoscopic radii), as well as 
direction-averaged one-dimensional radii are extracted. The radii decrease with kj, following a 
power-law behavior. This is qualitatively consistent with expectations from a collectively expanding 
system, produced in hydrodynamic calculations. The radii also scale linearly with (dVch/dr|)*^^. 
This behaviour is compared to world data on femtoscopic radii in heavy-ion collisions. While the de¬ 
pendence is qualitatively similar to results at smaller a decrease in the f?out/J^side ratio is seen, 

which is in qualitative agreement with a specific prediction from hydrodynamic models: a change 
from inside-out to outside-in freezeout conhguration. The results provide further evidence for the 
production of a collective, strongly coupled system in heavy-ion collisions at the LHC. 
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1 Introduction 


Collisions of lead ions at ^^^^=2.76 TeV have been recorded by A Large Ion Collider Experiment (AL¬ 
ICE) at the Large Hadron Collider (LHC) at CERN. In this energy regime, Quantum Chromodynamics 
(QCD) predicts the existence of a new state of strongly interacting matter, the quark-gluon plasma (QGP) 
in which quarks and gluons are no longer confined to individual nucleons. Experimental evidence for the 
existence of such matter has been found both at the Relativistic Heavy-Ion Collider (RHIC) |[ll-0] as well 
as at the LHC ll3-[iBl. The QGP behaves like a fluid with small viscosity and undergoes an explosive 
expansion. The study of the space-time structure and dynamics of this process may reveal important in¬ 
formation about the matter properties, such as its equation of state and the nature of the phase transition 
between deconfined and ordinary hadronic matter illLll2ll. 


Two-pion correlations at low relative momentum were first shown to be sensitive to the interaction vol¬ 
ume of the emitting source in p -|-p collisions by G. Goldhaber et al. 50 years ago [Il3n . Since then, they 


were studied in e+ -l-e lll4l] . hadron- and lepton-hadron il5l] . and heavy-ion 11614251] collisions. Espe 


cially in the heavy-ion case, two-particle femtoscopy has been developed into a precision tool to probe 
the dynamically-generated spatial structure of the emitting system. In particular, a sharp phase transition 
between the color-deconfined and confined states was excluded by the observation of short timescales. 
Moreover, femtoscopic measurement, together with other observations related to bulk collective flow, 
provide evidence that a strongly interacting system was created in the collision 112, 2^ 27]. 


Eemtoscopy in heavy-ion collisions is understood in some detail; for example see the experimental 
overview in 0] and model calculations in iEtUs^ . The dependence of the system size extracted from 
the data is investigated as a function of collision centrality and average transverse momentum of the pair 
kj = \pij + 7 ^ 2 ,T1/2. As the initial size of the system grows with increasing multiplicity (decreasing 
centrality), so does the apparent system size at freeze-out, measured by femtoscopy. Such increase is 
naturally produced in a hydrodynamic calculation. Strong hydrodynamic collective flow in the longitu¬ 
dinal and transverse directions results in the decrease of the apparent size of the system with increasing 
kj. This is because longitudinal and transverse velocity boosts cause particles emitted from spatially sep¬ 
arated parts of the collision region to move away from one another. Such particles cannot have a small 
momentum difference, and so correlation functions of boosted particles are sensitive to only part of the 
collision region. This part is referred to as the ’’homogeneity length” Olll . The decrease of the size with 


kj is observed in experimental data from heavy-ion collisions at all centralities, various collision energies 
and colliding system types, and is well described quantitatively in hydrodynamic models QS] and 


qualitatively in hadronic rescattering codes 1320. 


Taking into account the successful description of the femtoscopic scales at lower energies, the hydro- 
dynamic modelling has been extrapolated to collision energies of the LHC i30L 13314350. The expected 
increase in initial energy density (temperature) leads to larger evolution times, which in turn produce 
larger overall system size and stronger transverse and longitudinal flows. At the same time the freeze-out 
hypersurface evolves to have significant positive space-time correlation. This influences the radii of the 
system in the plane perpendicular to the beam axis. In particular the radius along the pair transverse 
momentum (called Rout) is decreased by the correlation with respect to the other transverse radius (called 
Rside). which decreases the Rout/Rside ratio All of those effects have been observed in the first mea¬ 
surement for central (0-5%) collisions at the LHC 1^ . This work extends this measurement to other 
centralities and compares the obtained radii to recent hydrodynamic calculations, in order to check their 
validity in a large range of event multiplicities. A measurement of one-dimensional radii was also per¬ 
formed using the two-pion and three-pion cumulants 071] . This work extends the two-pion measurement 
to several ranges of pair transverse momentum. 

The paper is organized as follows. In Section|2]the data taking conditions and data treatment is described. 
In Section[3]we give the details of the analysis of the correlation function. In Section|4|the extracted radii 
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are presented and compared to model expectations, while Section|5]summarizes our findings. 


2 Data taking and track reconstruction 


This work reports on the analysis of Pb-Pb collisions produced by the LHC during the 2010 datataking 
period. They were recorded by the ALICE experiment; the detailed description of the detector and the 
performance of all of its subsystems is given in |38,[^. Here we only briefl y d escribe the specific 


detectors used in this analysis. The ALICE Time Projection Chamber (TPC) ll40ll is a large volume 
gaseous ionization chamber detector, which was used both for tracking at mid-rapidity as well as for 
particle identification via the measurement of the specific ionization energy loss associated with each 
track. In addition to the TPC, the information from the ALICE Inner Tracker System (ITS) was used. The 
ITS consists of six cylindrical layers, two Silicon Pixel Detectors closest to the beam pipe, two Silicon 
Drift Detector layers in the middle, and two Silicon Strip Detectors on the outside. The information from 
ITS was used for tracking and primary particle selection, as well as for triggering. However the main 
triggering detector was the VO. It is a small angle detector consisting of two arrays of 32 scintillating 
counters. The first (VOA) is located 330 cm away from the vertex and covers 2.8 < r] < 5.1, the second 
(VOC) is fixed at the front of the hadronic absorber of the muon arm and covers —3.7 < r] < —1.7. The 
tracking detectors are located inside the solenoidal ALICE magnet, which provides a uniform magnetic 
field of 0.5 T along the beam direction. The TO detector 14111 was a main luminometer in the heavy-ion 
run. It consists of two arrays of 12 Cherenkov counters, covering —3.28 < fi — 2.97 and 4.61 < r] < 4.92. 
It has a time resolution of 40 ps. 


The minimum-bias trigger required a signal in both VO detectors, which was consistent with the col¬ 
lision occurring at the center of the ALICE apparatus. A total sample of approximately four million 
Pb-Pb events was used for this analysis. Lor the centrality range considered in this work the trigger effi¬ 
ciency was 100%. The centrality was determined by analyzing the signal from the VO detector with the 
procedure described in details in |[^. This ensured that the centrality determination was obtained using 
particles at significantly different rapidities than the ones used for the pion correlation analysis. This 
work presents results for seven centrality ranges: 0-5%, 5-10%, 10-20%, 20-30%, 30^0%, 40-50%, 
and 50-60% of the total hadronic cross section. The position of the event vertex in the beam direction 
with respect to the center of the ALICE apparatus was determined for each event. In order to ensure 
uniform pseudorapidity acceptance, only events with | 14 | < 8 cm were used in the analysis. 

Lor each event, a list of tracks identified as primary pions was created, separately for positive and negative 
particles. Each track was required to leave a signal both in the TPC and the ITS and the two parts of the 
track had to match. The TPC is divided by the central electrode into two halves, each of them is composed 
of 18 sectors (covering the full azimuthal angle) with 159 padrows placed radially in each sector. A track 
signal in the TPC consists of space points (clusters), each of which is reconstructed in one of the padrows. 
A track was required to be composed out of at least 80 such clusters. The parameters of the track are 
determined by performing a Kalman fit to a set of clusters. The quality of the fit is judged by the value 
of x^, which was required to be better than 4 per cluster (each cluster has two d.o.f.). The transverse 
momentum of each track was determined from its curvature in the uniform magnetic field. Two opposite 
field polarities were used through the data-taking period, for a check of systematic tracking effects. The 
momentum from this fit in the TPC was used in the analysis. Tracks which had a kink in the trajectory 
in the TPC were rejected. Trajectories closer than 3.2 cm in the longitudinal direction and 2.4 cm in 
the transverse direction to the primary vertex were selected to reduce the number of secondaries. The 
kinematic range for accepted particles was (0.14,2.0) GeV/c in transverse momentum and (—0.8,0.8) in 
pseudorapidity. Based on the specific ionization energy loss in the TPC gas dL/dx, a probability for each 
track to be a pion, kaon, proton, and electron was determined after comparing with the corresponding 
Bethe-Bloch curve. Particles for which the pion probability was the largest were used in this analysis. 
This resulted in an overall purity above 95%, with small contamination from electrons in the region 
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where the dE/Ax curves for the two particle types intersect. 

The accepted particles from each event are combined into same-charge pairs. The two-particle detector 
acceptance effects must be taken into account in this procedure. Two main effects are present: track 
splitting and track merging. Track splitting occurs when a single trajectory is mistakenly reconstructed 
as two tracks. The tracking algorithm has been specifically designed to suppress such cases. In a rare 
event when splitting happens, two tracks are reconstructed mostly from the same clusters in the TPC. 
Therefore, pairs which share more than 5% of clusters are removed from the sample. Together with the 
anti-merging cut described below this eliminates the influence of the split pairs. 


Two-particle correlated efficiency and separation power is affected by track merging. In the TPC, two 
tracks cannot be distinguished if their trajectories stay close to each other through a significant part of 
the TPC volume. Although this happens rarely, such pairs by definition have low relative momentum 
and therefore their absence distorts the correlation function in the signal region. The effect of track 
merging has been studied in central collisions in the previous work Ool . In this work we have used 
a similar procedure to correct for the merging effects, through dedicated two-particle selection criteria. 
More details are given in Sec. 13.31 


3 Correlation function analysis 

The two-particle distribution for same-event pion pairs depends on several factors, including trivial 
single-particle acceptance effects. To extract only the relevant two-particle correlation effects, the corre¬ 
lation function formalism, described below, is applied. 

3.1 Correlation function construction 

The femtoscopic correlation function C is constructed experimentally as 


C{q) = 


Mg) 

B{qy 


( 1 ) 


where q = p\ — p 2 is the pair relative momentum (due to fixed masses of the particles only three com¬ 
ponents are independent). The magnitude of this vector is referred to as ^inv. For a detailed description 
of the formalism, see e.g. 1421] . The signal distribution A is composed of pairs of particles where both 
come from the same event. The background distribution B is constructed with the “mixing” method, in 
which the two particles come from different events, which must have similar characteristics, so that their 
single-particle efficiency and distribution are as close as possible. To form a “mixed” pair, particles must 
come from two events, for which the centralities differ by no more than 2.5% and vertex positions differ 
by no more than 4 cm. The correlation function is normalized with the ratio of the number of pairs in 
the B and A samples in the full q range used (0 — 0.25 GeV/c), so that C at unity means no correlation. 
The dependence on pair momentum sum is investigated by doing the analysis for various ranges of kj, 
namely: 0.2-0.3, 0.3-0.4, 0.4-0.5, 0.5-0.6, 0.6-0.7, 0.7-0.8, and 0.8-1.0 GeV/c. The ranges were the 
same for each centrality range, which overall resulted in 49 independent correlation functions per pair 
charge combination. 


The momentum difference q is calculated in the Longitudinally Co-Moving System (LCMS), in which 
the pair total longitudinal momentum vanishes: pi,L + P 2 ,L = 0. The three coordinates of q in LCMS are 
defined as follows: long - along the beam axis, out - along the pair transverse momentum, and side - 
perpendicular to the other two. In Fig.[T]the projections of three example correlation functions along these 
axes are shown. A significant, approximately Gaussian enhancement at low relative momentum is seen 
in all projections. The width of the correlation grows with increasing cenhality (lowering multiplicity) 
as well as with increasing ^t. 
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Fig. 1: Projections of the Cartesian representation of the three-dimensional two-pion correlation functions along 
the out (a), side (b) and long (c) axes. The centrality and pair momentum ranges for the three functions are given 
on the plot. In each case the other components are projected over ± 20 MeV/c around 0 in the other q directions for 
central collisions at low kj and ± 33 MeV/c in the other two cases. Lines represent the corresponding projections 
of the ht to the experimental correlation functions with formula from Eq. (Ell. 


The pair distributions for identical particles have specific symmetries, which are naturally represented 
in a spherical harmonic (SH) decomposition 1431 Udll. In particular for identical particles all odd-Z and 
odd-m components of a correlation function representation vanish. The gross features of the correlation 
function which are relevant for femtoscopy are fully captured by the low-Z components of the decompo¬ 
sition. The I = 0,m = 0 component is sensitive to the overall size of the pion source. The I = 2,m = 0 
component is sensitive to the difference between longitudinal and transverse size, and the I = 2,m = 2 
component reflects the difference between sideward and outward component of the transverse radius. 
Therefore, those three components of the SH representation contain the same information as the Carte¬ 
sian one for the purpose of the femtoscopic analysis. In particular both representations allow for fitting 
of the correlation function with the same theoretical formula. The next non-vanishing components are 
for Z = 4. Their analysis is beyond the scope of this paper, which focuses on the overall width (variance) 
of the distribution in three directions. 


In Fig. |2] we show the first three non-vanishing components of the spherical harmonics representation 
of three example correlation functions, the same as in Fig. [T] In the (0,0) component the enhancement 
at low-^ is clearly visible and its width is increasing with centrality and kj. In this representation the 
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Fig. 2: First three non-vanishing components of the SH representation of the two-pion correlation functions, 
I = 0,m = 0 in (a), I = 2,m = 0 in (b) and / = 2,m = 2 in (c). The centrality and pair momentum ranges for the 
three functions are given on the plot. The lines show the corresponding SH components of the fit with formula 
from Eq. dUl. 


low-^ dip coming from Coulomb repulsion is visible better than in the Cartesian one. For large Ut;, 
there is a lack of low-i^ pairs, as a result of track merging (see below). Nevertheless, the correlation 
clearly extends well beyond the region with low statistics. The other two components: (2,0) and (2,2), 
show a non-trivial correlation structure deviating from zero, indicating that the shape is not spherical, but 
rather ellipsoidal in LCMS. The width of these structures is consistent with the width of the correlation 
enhancement in the (0,0) component, which means that all three reflect the properties of the femtoscopic 
signal. The analysis of the shape of these structures is the main focus of the next section. 

If the available statistics is limited, as is sometimes the case for lower collision energies or particles 
heavier than pions, the analysis is performed only as a function of magnitude of relative momentum 
most naturally calculated in the Pair Rest Frame (PRF). In this work we present results in this variable 
for completeness. Data was analyzed in the same centrality and pair k-j ranges as the ones used in the 
three-dimensional analysis. 
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3.2 Fitting the correlation function 

The space-time characteristics of the source are reflected in the correlation function. They are connected 
via the Koonin-Pratt equation 

Ciq)= Js{r,q)\^{lr)\^d\ (2) 

where r is the pair space-time separation four-vector. S is the source emission function, interpreted as a 
probability to emit a pair of particles with a given relative momentum and space-time separation. *P is 
the two-particle interaction kernel. In the simplest case of non-interacting particles (e.g. photons) it is 
the modulus of the pair wave-function. If the Coulomb or strong interaction between the particles (called 
Final State Interaction - ESI) needs to be taken into account, then 'P becomes the Bethe-Salpeter am¬ 
plitude corresponding to the solution of the relevant quantum scattering problem, taken with the inverse 
time direction ll45n . 

Previous studies at RHIC 116-18, 22, 23, 4^ 471 and at the EHC have approximated the source by 
a Gaussian, treating any difference between the real data and a Gaussian as a correction. This procedure 
was also universally used in all past pion femtoscopic analyses of heavy-ion collisions. Therefore, we 
also use it here by writing 


S{r,q) 


exp 


' out 
4/?2 


r2. 

side 

4/?2 

^^side 


long 


(3) 


where rout, Aide, and riong are components of the relative separation r. This static form of S is expressed 
in ECMS, with Rout, ^side, and Riong being the single-particle source sizes of the system later referred to 
as “femtoscopic radii”, or simply “radii”. They quantify the lengths of homogeneity of the system in the 
outwards, sidewards, and longitudinal direction, respectively. 


For like-sign pions the strong interaction contribution is small for the source sizes expected here (a 
few fm) 142(1, so it is neglected. The remaining 'P is a convolution of the Coulomb interaction and 
wave-function symmetrization. As an approximation, the Coulomb part is factorized out and integrated 
separately in the procedure known as the Bowler-Sinyukov fitting i48(l49ll . It is well tested and is appli¬ 
cable for pions and for large source sizes expected in this analysis. In this approximation the integration 
of Eq. ^ with S given by Eq. Q gives the following fit form for the correlation function 


C{q) = N{\-X) (4) 

-|- AA,R(;;(^inv) [l “hCXp ( ^out^out -^side^side ^long^long)] ) 

where N is the overall normalization factor. The function Rc(<?inv) is the Coulomb part of the two-pion 
wave-function integrated over the spherical Gaussian source with a given radius. For each correlation 
function it is set to the value from the one-dimensional analysis (see Section 1431) . to reflect the decrease 
of the source size with multiplicity and ^t. Its variation is a source of systematic uncertainty. The dilution 
parameter X is introduced to account for the fact that not all measured pion pairs are correlated, and that 
the real emission function may deviate from a Gaussian form. 

The fit is performed with the log-likelihood method for the three dimensional correlation function in 
Cartesian representation, resulting usually in several thousand degrees of freedom. Examples are shown 
in Fig. [U The Gaussian fit is able to reproduce the overall width of the correlation in all cases. Some 
details of the behavior at low-^ may not be perfectly described, which can be attributed to the limitations 
of the Bowler-Sinyukov formula as well as to the non-Gaussian, long-range tails which may be present in 
the source. Some deviations from the Gaussian ellipsoid shape for the higher centrality can also be seen 
for the long direction. We leave the detailed investigation of these effects for future work. Nevertheless, 
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Uncertainly source 

Rout [%] 

^side [%] 

■^long [%] 

CE represenlafion 

0.5-5 

0.5-4 

0.5-8 

Dalasel comparison 

< 1.5 

< 1.5 

<2 

Eil-range dependence 

0.5-4 

0.5-3 

1-5 

Two-frack cul variafion 

3-10 

2-12 

2-13 

Coulomb correcfion 

3 

1 

1 

Momenfum resolution correcfion 

2 

2 

2 

Cenfralify esfimalion 

1 

1 

1 

Tolal 

6-14 

4-13 

4-17 


Table 1: List of contributions to the systematic uncertainty of the extracted femtoscopic radii (ranges given). 


the overall sizes R of the system, which are mostly sensitive to the width of the correlation are well 
estimated. The deviation of the correlation function from the pure Gaussian shape is smaller than a 
similar deviation in pp collisions tSfll . 


An equivalent fit is also performed for the SH representation of the correlation. Eq. ([Hi is numerically 
integrated on a (p — 9 sphere for each ^lcms bin, with proper F/” weights, to produce the 3 components 
of the SH decomposition. Statistical uncertainties on each component are taken into account, as well 
as the covariance matrix between components. Examples are shown in Eig. |2] The fit describes the 
general direction-averaged width of the correlation function, shown in the upper panel. Small deviations 
can be seen in the shape and the behavior of the fit at small q, as discussed earlier for the Cartesian fit. 
The deviations from zero in the (2,0) and (2,2) components are small but statistically significant for 
a large number of multiplicity/kx ranges, indicating that the source size is slightly different in all three 
directions. 


Eor the one-dimensional correlation functions constructed as a function of ^inv, the fit was performed 
with a simplified version of Eq. ([H) 


C{q) = {\-X) (5) 

+ X/fc [1 + 

where fhe only radius paramefer is /?inv - the one-dimensional direction-averaged femtoscopic radius in 
the PRE. 

3.3 Systematic uncertainties of the radii 

Tabled] lists the systematic uncertainty contributions. The range of values is given to provide a general 
estimate of the importance of each contribution, however the systematic uncertainty is estimated for each 
point individually. Separate analyses were performed for positive and negative pions, as well as for two 
datasets collected with opposite polarities of the magnetic field inside AEICE. This resulfs in four inde- 
pendenf dafa samples for which cerfain sysfemafic effecls, mosf nofably fhe single-frack inefficiencies, 
are differenf. Correlafion functions for all four samples for all cenfralify and pair momenfum ranges are 
sfalislically consisfenf, after all corrections are applied; fhis is an imporfanf systematic cross-check of 
fhe mefhodology. In fhe following discussion fhe cenfral values are sfafislical averages of fhe til values 
obfained for fhe four samples. The sysfemafic uncerfainfy arising from differences among fhe dafa sam¬ 
ple is befween 1 and 2% for all radii. The ofher sysfemafic uncerfainfies are analyzed for each sample 
separately; fheir final value is fhe convolufion of fhe uncerfainfies for each sample. 

Two correlafion funcfion represenfafions are used in fhis work: fhe Carfesian and spherical harmon¬ 
ics. They are mafhemafically equivalenf, fhe tiffing procedure used fhe same functional form for bofh. 
However fhe implemenfafion of fhe tiffing procedure is quife differenf: log-likelihood vs. regular 
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fit, three-dimensional Cartesian histogram vs. three one-dimensional histograms, fitting range as three- 
dimensional cube in <?out)^side)<?iong or a three-dimensional sphere with constant ^lcms radius among 
others. Therefore, the fits to the two representations differ systematically upon variation of the fitting 
procedure (fit ranges, Bowler-Sinyukov approximation, etc.). The difference between the values for the 
two fits is taken as a part of the systematic uncertainty. It usually ranges from 1% to 3% and grows with 
pair k-j and multiplicity. 

Variation of single particle cuts around the default value results in modifications of single particle ac¬ 
ceptances and purities. However the correlation function shape should be to the first order insensitive to 
those effects. We have checked that for a reasonable variation of the single particle cuts, the resulting 
radii are consistent within statistical uncertainties. 


The measurement of the average event multiplicity for a given centrality range has a known uncertainty 
of 3-4% for all centrality classes fl]. The femtoscopic radii in heavy-ion collsions were observed to 
scale linearly with (dVch/dri)^'^^at lower collision energies |20]. Such trend is also predicted by hydro- 
dynamical models and is expected to hold at LHC. Therefore the systematic uncertainty coming from 
the multiplicity estimation is about 1%. 

The dominant systematic effect on the two-particle correlation function is the two-track correlated effi¬ 
ciency. This effecf has been sfudied in previous work 061] for cenfral collisions. The effecl of “spliffing” 
is fully removed wifh fhe help of dedicated fwo-frack selecfion criferia, mentioned in Sec. |2j and does 
nol influence fhe radii. However when fwo particles’ frajecfories are locafed close fo each ofher in fhe 
volume of fhe TPC defecfor, fhey may be reconsfrucfed as one frack or nol be reconslrucled al all. This 
effecl, called “merging” in fhe following, resulls in a loss of reconslruclion efficiency for such pairs of 
fracks. Pairs of primary pion frajecfories close in space correspond lo low relalive momenlum, fhere- 
fore “merging” will affecl fhe pion correlalion function exaclly in fhe femtoscopic signal region. The 
Iwo-lrack efficiency was sfudied in Monfe-Carlo simulation of fhe ALICE defecfor. Lor pion pairs al low 
relative momenlum an efficiency loss of up to 20% was observed. A Iwo-lrack selection was chosen, 
such lhal fhe resulting correlalion function was nol affecled by fhe inefficiency. “Merging” affecls only 
pairs which are spatially close in fhe defector. The “closeness” can be quanfified by fhe pseudorapidify 
difference Ar] for fhe pair and only pairs wifh |Ari| < 0.016 are affecled. The frajecfories musl also be 
close in fhe Iransverse plane, where fhey are curved by fhe magnelic field. The azimulhal coordinate tp* 
of fracks al a radius of 1.2 m from fhe collision poinl (so roughly in fhe center of fhe TPC volume) is 
calculated. Pairs wifh Y^Ap^ + Atp^ < 0.045 are affecled by merging. Pairs which simullaneously sal- 
isfy Ihis and fhe previous angular criteria are removed bofh from fhe signal and fhe background sample. 
The correlation funclion calculated for pairs surviving Ihis cul is Ihen nol affected by “merging”. As a 
sysfemafic check a differenl procedure fo calculale tp* was used, where fhe value was laken nol al a fixed 
radius, buf inslead al a radius inside fhe TPC (i.e. belween 0.5 and 2.4 m) where (p* was fhe smallesl. 
The differences belween radii for fhe correlation functions calculated wifh fhe fwo mefhods is laken as 
sysfemafic uncerfainfy. If is from 2-13%. 

A pair of same charge pions Iraversing a solenoidal field can have fwo configuralions: a “sailor” where 
the two trajectories bend away from one another and never cross, or a “cowboy” where the trajectories 
bend toward each other and can cross lllSn . Merging affects the “cowboy” pairs, but only weakly influ¬ 
ences fhe “sailors”. The fwo configuralions correspond to differenl phase space regions in fhe Cartesian 
represenlafion of fhe correlalion function. We have performed fwo fils fo fhe correlations (correcled for 
merging), resfricfing fhe filling range lo eilher fhe “cowboy” or fhe “sailor” region. The consistency 
of fhe radii obfained in Ihese fils was used as an esfimafor of fhe effecliveness of fhe cul for removing 
merged fracks. A less restrictive cut degraded agreement between the “cowboy” and “sailor” radii, while 
making the cut stricter did not improve the agreement but reduced statistics in the signal region. The 
comparison allowed optimization of the cut, and provided another estimate of the systematic uncertainty 
on the two-track correction procedures. The uncertainty estimated in this way is consistent with the un- 
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certainty determined by varying the tp* definition. This uncertainty is largest for transverse radii, is most 
prominent for high multiplicity events (central collisions), and affects a wider region in q for pairs with 
higher 

The separation into the “cowboy” and “sailor” phase-space regions is not feasible for the SH represen¬ 
tation of the correlation function. In this case, when the anti-merging cut is not properly applied, one 
observes significant non-femtoscopic signals, especially in the C\ component of the correlation function 
and a reliable fit cannot be performed. Therefore, the Cf component serves as a sensitive independent 
check of the effectiveness of the “anti-merging” cut. This is a good illustration on how the two represen¬ 
tations complement each other in the systematic studies. 

The fit was performed for several values of the fitting range in q (varying with multiplicity and pair kj, 
following the changes in the correlation function width). The variation of the fitted radii with the change 
of the range was taken as another component of the uncertainty, which is less than 5% for all the radii. 


In addition to the uncertainties listed above, other systematic effects can influence the extracted radii. 
The first is the momentum resolution, which was studied in iB^ . The correction procedure described 
there is used in this work, as well. The uncertainty on the radii from this correction is 2%. Another 
effect is the influence of the Bowler-Sinyukov procedure on the extracted radii and X parameter (fraction 
of correlated pairs). The procedure results in an uncertainty of 3% on /?out^ 1% on the other radii, and 
lowers the X value by up to 5% 12811 . 


All the systematic uncertainty components mentioned above are added in quadrature, the range of values 
of the total systematic uncertainty is given in Table [T] 


The one-dimensional analysis is performed in PRF, where the total momentum of the pair vanishes. In 
the transformation from LCMS to PRF Rout is scaled by the y factor for the pair, depending on kx and 
as a consequence is larger than the other two components. In such case the direction-averaged one¬ 
dimensional correlation function becomes non-Gaussian. This produces a dependence of the fit value on 
the range of the fit, resulting in a systematic uncertainty of up to 10%. Other components of the uncer¬ 
tainty, such as field orientation dependence, momentum resolution and Coulomb correction dependence 
are comparable to those from the three-dimensional fit. 


4 Results 

4.1 Three-dimensional radii 

The outcome of the fitting procedure are 49 sets of femtoscopic radii, one set for each centrality and 
kj range. They are shown in panels a)-c) of Fig. [3] The radii in all directions are in the range of 2 to 8.5 
fm. The radii universally decrease with increasing kx, in qualitative agreement with a decreasing homo¬ 
geneity length, as predicted by hydrodynamics. Such behavior is a strong indication of a large degree of 
collectivity in the created system. The radii are also universally higher for more central collisions, which 
correspond to growing final-state event multiplicity. For the lowest kx, Riong is generally the largest, 
whereas at large kx there is no universal ordering of the radii. 

In the most central collisions the values of the X parameters are around 0.33 for the lowest kx range, are 
increasing linearly to 0.43 for k-j of 0.65 GeV/c, and falling to 0.36 in the last k-j range. For peripheral 
collisions, they are higher by 0.05 to 0.07. The lowering of X from a theoretical maximum of 1.0 can be 
attributed to a number of factors. The pion sample contains daughters of long-lived strongly decaying 
resonances; their fraction falls with growing px- Asa result the source function contains non-Gaussian 
tails extending to large relative separation. There may be additional factors influencing the shape of the 
correlation function, coming from the source dynamics. The exact source shape usually deviates from a 
Gaussian, in a way that lowers the X parameter of the Gaussian fit significantly. The detailed study of this 
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Fig. 3: Femtoscopic radii as a function of pair transverse momentum kj for seven centrality ranges. /?out is shown 
in (a), Rside in (b), f^iong in (c) and /?out/f^side in (d). The points for centralities 5-10%, 10-20%, 30^0% and 
40-50% have been slightly shifted in the x direction for visibility. Statistical uncertainties are shown as black lines 
(usually smaller than symbol size), systematic uncertainties are shown as colored caps. 


shape requires a dedicated methodology and is beyond the scope of this work. In addition, at the largest 
kj the electron and pion dE /dr become comparable in the TPC, and the pion sample is contaminated by 
electrons, which also lowers X. The approximate treatment of the Coulomb interaction in the Bowler- 
Sinyukov fitting procedure lowers it by 5-10% ll28n . This effects is most pronounced for large source 
sizes (central collisions). Finally, possible coherent emission of pions 115IH is expected to lower X by a 
few percent. The X parameter values given here are for reference only and their physics interpretation is 
not discussed. 


In panel (d) of Fig.|3]the /?out/^side ratio is shown. Its systematic uncertainty is determined independently 
from those of R^ut and /?side^ to account for the fact that they may be correlated. The ratio is consistent 
with unity for central collisions. Its value slowly decreases for more peripheral collisions, and reaches 
0.85 for peripheral collisions and high A:t. Based on hydrodynamic models, the f?out/^side ratio has 
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been proposed as a sensitive probe of the shape and space-time correlation present at the freeze-out 
hypersurface OSL15211 . In particular this ratio at the LHC was predicted to be lower than a value of 1.1 
measured at top RHIC collision energies. 
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Fig. 4: The radii plotted as a function of (chVch/dp)^^^ for seven kj ranges. f?out is plotted in (a), f?side in (b), and 
/?iong in (c). The dashed lines represent linear fits to each set of points for a given component and kj range. Points 
for some kj ranges have been slightly shifted in the x direction for visibility. Statistical uncertainties are shown as 
black lines (usually smaller than symbol size), systematic uncertainties are shown as colored caps. 


4.2 Scaling of the radii 

It has been argued in |[^ that the femtoscopic volume scales with the final-state event multiplicity, and 
that each of the three-dimensional radii separately scales with this value taken to the power 1 /3. In Fig.|4] 
we present the dependence of the radii on multiplicity for Pb-Pb collisions. The scaling is evident for all 
datasets, for all three directions and all analyzed pair momentum ranges. 


Similarily, hydrodynamics predicts approximate scaling of the radii with pair transverse mass mx = 
+ 135]. The slope parameters of the lines shown in Fig. |4]are plotted in Fig. |5l as a function of 
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Fig. 5: The slope parameters of the linear fits shown in Fig.|4l as a function of pair mj. Lines represent the power- 
law fits (see text for details). Errorbars represent the uncertainty of the parameter of the hts, in which the combined 
systematic and statistical uncertainties on the radii ware taken into account. 


mj. They are fitted with a power law function of the form 


( 6 ) 

\m%J 

where P and a are free parameters. The slope parameters follow the power law scaling within the current 
systematic uncertainties, the value of the a parameter is —0.65 zh0.12 for long, —0.46 ±0.13 for out, 
and —0.52 ±0.11 for side direction. The dependence of the values of femtoscopic radii on centrality and 
ki factorizes into a linear dependence on (dAch/dr])^^^ and a power law dependence on ott. 

4.3 One-dimensional analysis 
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Fig. 6: One-dimensional femtoscopic radius R^ny as a function of pair transverse momentum kj for seven centrality 
ranges. Points for 5-10%, 20-30%, and 40-50% centrality have been slightly shifted in the x direction for visibility. 
Statistical uncertainty is shown as black lines, systematic uncertainty is shown as colored caps. 
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The results of the one-dimensional fits with Eq. ® are shown in Eig.[^ Similarly to the three-dimensional 
case, the radius is increasing with event multiplicity (decreasing centrality). Therefore, the final state 
shape is reflecting the growth of the initial shape with decreasing centrality. The is also decreasing 
with pair transverse momentum. This is usually understood as a manifestation of the hydrodynamic col¬ 
lectivity. The one-dimensional radius also serves as a comparison basis with the femtoscopic analysis 
for heavier particles, where the one-dimensional analysis is standard and the three-dimensional analysis 
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* STAR Au-Au = 200 GeV 

O STAR Cu-Cu = 200 GeV 

T STAR Au-Au = 62 GeV 
0 STAR Cu-Cu = 62 GeV 

* CERES Pb-Au ^ = 17.2 GeV 

* ALICE Pb-Pb = 2760 GeV 


Fig. 7: Comparison of femtoscopic radii, as a function of measured charged particle multiplicity, for a number 
of collision systems and collision energies IJJ, 4^ 53]. Dashed lines show linear fits done to heavy-ion data, 
excluding ALICE (dotted lines show one sigma contour). Various experiments use differing ky ranges, values 
shown on the plot are for the range for which the average ky is closest to the selected value of 0.4 GeV/c; in case of 
ALICE an average of two neighbouring ranges of k-y is shown. Systematic uncertainties are shown where available. 


is challenged ^ the more complicated description of the pair interaction as well as significantly smaller 
statistics ll54l55l]. 


4.4 Comparison to previous measurements 

In Fig. |7] the heavy-ion data from Pb-Pb collisions at the LHC reported in this work are compared to the 
previous measurements, including results obtained at lower collision energies. It has been argued 1201] 
that three-dimensional femtoscopic radii scale with the cube root of measured charged particle multi¬ 
plicity not only for a single energy and collision system, but universally, across all collision energies and 
initial system sizes. The dashed lines in the figure are linear fits to heavy-ion data available before the 
startup of the LHC (the dotted lines show one sigma contours of these fits). At lower energies the linear 
scaling was followed well in long and side directions and only approximately in out, with some outliers 
such as the most peripheral collisions reported by STAR. Our data at higher collision energy show that 
the scaling in the long direction is preserved. The data for the side direction fall below the scaling trend, 
although still within the statistical uncertainty. A clear departure from the linear scaling is seen in the 
out direction; data from the LHC lie cl early below the trend from lower energies. Such behavior was 
predicted by hydrodynamic calculations 0311 and is the result of the modification of the freeze-out shape. 
Larger initial deposited energy produces larger temperature gradients and longer evolution time at LHC. 
This results in a change from outside-in to inside-out freeze-out and this modification of the space-time 
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correlation drives the /?out/^side ratio to values lower than at RHIC. Therefore, already for heavy-ion data 
the simple universal linear scaling is broken in the transverse direction. As already reported in 15011 . the 
femtoscopic radii for pp collisions also exhibit linear scaling in the same variable, albeit with signifi¬ 
cantly different parameters. In this case the scaling between different colliding systems is broken again 
in the longitudinal direction. 

Other scaling variables were proposed for the femtoscopic radii, based on Monte-Carlo simulations of the 
initial state, such as average number of participants Apart 14611 or the characteristic initial transverse size 
R 1^ 57]. If data from different centralities but the same collision energy are considered, a linear scaling 
is indeed observed for radii in all directions, if they are plotted as a function of either of these variables. 
However, no such scaling is observed when data from two collision energies (e.g. top RHIC and top 
LHC) is compared. In that sense these variables are less adequate than (dAch/dr]) for which the linear 
scaling is preserved, across many collision energies, for at least one direction (Riong)- This observation is 
consistent with an expectation that the final freeze-oul volume, reflecled in fhe femtoscopic radii, should 
scale wifh fhe final-sfafe observable (such as e.g. (dAch/dr])^^^), while fhe simple “geomefric” initial 
sfafe variables do nol confain enough information. They musf be replaced wifh modeling of fhe full 
evolution process, such as fhe one given by hydrodynamics, which depends on addifional paramefers 
apart from fhe inifial size, such as initial energy densify. 


4.5 Model comparisons 


The hydrodynamic models predicf bofh fhe cenfralify and pair momenfum dependence of fhe femtoscopic 
radii. Usually fhe paramefers of fhe model (inifial energy densify profile, fhe equafion of sfafe and 
fhe freeze-ouf condifion) are adjusfed fo reproduce fhe shape of fhe single particle inclusive fransverse 
momenfum specfra as well as fhe elliptic flow. The charged parficle mulfiplicify musf also be reproduced, 
preferably as a function of pseudorapidify if fhe model employs full fhree-dimensional modeling. The 
fofal particle mulfiplicify is usually defermined on fhe freeze-ouf hypersurface by employing sfatislical 
hadronizafion. Affer converting fhe confinuous medium fo hadrons, final sfafe inferacfions are faken info 
accounf eifher in fhe simplified form, wifh propagafion and decay of hadronic resonances or wifh fhe full 
rescaffering simulafion. We compare our resulfs to predictions from fhe Therminafor model coupled fo 
(3-1-1)D viscous hydrodynamics 134, 35, 58]. Similar resulfs have been obfained in fhe HKM model 0]. 


In Fig. [H we show fhe comparison of our dafa to fhe calculafions from fhe (3-1-1 )D hydrodynamic model 
coupled to fhe Therminafor sfafisfical hadronizafion code. The model is fully fhree-dimensional and if is 
able fo reproduce values of Riong for all cenfralifies, wifh some overpredicfion of fhe overall magnifude 
and fhe slope of fhe kj dependence, especially af low momenfum. This is an indicafion fhaf fhe longi- 
fudinal dynamics is reasonably described in fhe model, bofh in momenfum and space-time sectors. The 
Rout is well described for all cenfralifies, bofh in magnifude as well in fhe slope of fhe kj dependence. 
The slope of fhe kj dependence is also well described for Rside> buf fhe magnifude is lower fhan in dafa, 
allhough wifhin fhe syslemalic uncertainly. The infercepf of Rside af low kj is usually associafed wifh 
fhe overall geomelrical size of fhe syslem, while fhe slope of fhe kj dependence of bofh fransverse radii 
depends on fhe amounl of flow in fhe syslem. Bofh are well reproduced, so fhe hydrodynamic approach 
is in good agreemenl wifh our dafa. 

Anolher model based on hydrodynamic formalism, fhe HKM is also shown in Fig. [8] for 0-5% 
mosl cenlral collisions. If differs from fhe previous in fhe implemenfafion of fhe freeze-ouf process. If 
also direclly Ireals hadron rescatfering wifh fhe UrQMD simulafion. Neverlheless, fhe pion femloscopy 
in cenlral Pb-Pb collisions al fhe LHC is reproduced in fhe calculation, wifh some underprediclion of 
Rside and Riong magnifude. Therefore, fhe approximafe agreemenl wifh dafa is a universal fealure of 
such models, nof of a particular implemenfafion. The particular choice of inifial condilions and fhe 
equation of sfafe for Ihese models was molivaled by fhe analysis of RHIC femtoscopic dafa. This choice is 
essenfial for fhe correcl descripfion of fhe dafa, due to fhe collecfive flow mechanism in fhe hydrodynamic 
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Fig. 8: Comparison of the femtoscopic radii (Rout in (a), Rside in (b), and /?iong in (c)), as a function of pair trans¬ 
verse momentum, with the calculation from the Therminator and (3 h- 1)D hydro model 1^ . for four centralities 
(here identihed by the value of the impact parameter b used in the calculation) and with the HKM model for 
the central data. Closed symbols are experimental data (with statistical and systematic uncertainty), bands and 
dashed lines are the model calculations. 


evolution at the LHC energies. It also indicates that the details of the freeze-out process have limited 
influence on pion femtoscopy. Some studies suggest that femtoscopy of heavier particles might be a 
more sensitive probe in this case BOT . 


The data to model comparison of the f?out/^side ratio is plotted in Fig. |9l It shows values consistent 
with unity for central collisions, both for models and data. Such low values are associated with the 
change to outside-in freeze-out scenario at LHC collision energies. For most peripheral collisions the 
ratio decreases with kj even more, to values smaller than unity. This decrease is qualitatively reproduced 
in the model, although the calculations are at the upper edge of the experimental systematic uncertainty. 
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Fig. 9: Comparison of the /?out/f^side ratio, as a function of pair transverse momentum, with the calculation from 
the Therminator and (3 h- 1)D hydro model 03511 . for 0-5% centrality in panel (a) and 50-60% centtality in panel 
(b). The comparison to the HKM model S is also shown for the central data. Closed symbols are experimental 
data (with statistical and systematic uncertainty), bands and dashed lines are the model calculations. 


5 Conclusions 

We reported on the centrality and pair dependence of the three-dimensional and direction-averaged 
one-dimensional pion femtoscopic radii in Pb-Pb collisions at y^5NN=2.76 TeV. The behavior of the 
femtoscopic radii can be factorized into a linear dependence on cube root of charged particle density 
(separately for each of the three-dimensional radii) and a power-law dependence on pair transverse mo¬ 
mentum, with slightly different exponents for each direction. The dependence for the longitudinal radius 
is steeper than for the transverse ones. The dependence was also compared to the one observed in 
heavy-ion collisions at lower energies and to other collisions systems. The radii at the LHC follow the 
“universal” (dAch/drj)'^^ scaling in the long direction, but the radii in the transverse directions are 
below the “universal” curve. Simple linear scaling predictions are not valid when the collision energy 
is increased by an order of magnitude. The details of the dynamic evolution of the system influence 
the results significantly. This is in qualitative agreement with predictions from hydrodynamic models. 
In particular, when moving from RHIC to LHC collision energies, they produce a change in freeze-out 
shape, larger transverse radial flow and longer system evolution time. Comparison of the full dataset 
to the calculations from the recent hydrodynamic models, including three-dimensional evolution as well 
as hadronic stage, generally show a good agreement, which is complementary to similar agreement ob¬ 
served for momentum-only observables, such as momentum spectra and elliptic flow. The existence 
of such agreement both in the space-time as well as in momentum sectors, provides strong arguments 
for the validity of hydrodynamic models for the description of flowing bulk matter created in heavy-ion 
collisions at the LHC. 
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